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ABSTRACT 



We study the thermodynamics of massive Gross-Neveu models with ex- 
plicitly broken discrete or continuous chiral symmetries for finite temperature 
and fermion densities. The large N limit is discussed bearing attention to the 
no-go theorems for symmetry breaking in two dimensions which apply to the 
massless cases. The main purpose of the study is to serve as analytical ori- 
entation for the more complex problem of chiral transition in 4— dimensional 
QCD with quarks. For any non-vanishing fermion mass we find, at finite 
densities, lines of first order phase transitions. For small mass values traces 
of would-be second order transitions and a tricritical point are recognizable. 
We study the thermodynamics of these models, and in the model with broken 
continuous chiral symmetry we examine the properties of the pion like state. 



1. Introduction 



We have examined the massive Gross-Neveu model |IJ at finite temperature and density 
within the mean field approximation. In spite of the problems related to the low dimen- 
sionality of the model, this may represent in our opinion a guide to the thermodynamics 
of chiral symmetry restoration in QCD. We have studied the model with a bare mass 
term included, always kept non vanishing. Due to such a choice, no chiral phase transi- 
tion is present, since chiral symmetry is explicitly broken from the beginning and thus, 
strictly speaking, the Mermin- Wagner- Coleman theorem H@ does not apply. Although 
the symmetry is explicitly broken we find, as in our previous study for a QCD model 
H, that some first order phase transition still survives. The critical line obviously moves 
in the plane of temperature and chemical potential for growing masses, but it survives 
even for large mass values. We derive the equation of state and study the phase dia- 
gram for different choices of thermodynamical variables. We discuss the isotherms in 
the pressure- inverse density plane, which resemble the Van der Waals isotherms for the 
vapour-liquid transition in the water. The construction is made through the study of 
the effective potential, which contains all the physical information about the stable and 
metastable phases of a given model. We can define a critical point as the ending point 
of the coexistence region. There are, below this point, two regions of very low and of 
very high compressibility, separated by the coexistence region which ends at the critical 
point. 

We have also considered the explicit (small) breaking of a continuous symmetry in or- 
der to implement the study of soft-pion type properties in the model. It is evident that 
(differently from 4— dimensional QCD) in two dimensions the pion-like particle cannot 
be considered as a Goldstone particle due to the Coleman theorem. Thus one principal 
difference between previous studies [[| and the present work is that here the zero mass 
limit cannot be taken without a complete change of the picture. 

To further comment on how no-go theorems work in a finite temperature field theory and 
to better motivate our presentation, we start by summarizing in a very schematic way 
known results in Section 2. 

Section 3 is devoted to the results for the phase diagram and the equation of state of the 
model with a broken discrete symmetry. Here also naive zero mass limit results will be 
presented, to better clarify those for the massive case in which we are indeed interested, 
and at the same time to show what would be reasonable to expect in analogous studies 
in 3+1 dimensions. In Section 4 we attempt at a description of pion properties in the 
model with a broken continuous symmetry. As the effective potential in the mean field 
approximation can be put in the same form as that of Section 3, the results concerning 
the equation of state and the phase diagram are the same. Finally, some useful calcula- 
tions are summarized in the Appendix. 
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2. General review 



The Gross- Neveu model [|IJ is a well-known two-dimensional theory with four-fermion in- 
teractions which is asymptotically free. The fermion field has iV components. The model 
was originally considered in the 1/N expansion. The massless formulations, with discrete 
or continuous symmetry, have been extensively studied for zero or for finite temperatures 



and densities |p|-||12||, giving rise to several discussions related to the low dimensionality 
of the model, with its implications for symmetry breaking and phase transitions, and to 
the validity of the 1/N expansion (many aspects have been already discussed in ref. (T! 



Discrete symmetry 

Let us summarize in a very schematic way the main known results. The Lagrangian is 

C = ^ + ^g 2 (^) 2 (2.1) 

(for studies of the limit N —>■ oo one will also define g 2 N = A). It is invariant under the 
discrete chiral symmetry 

i> -> 75^ (2.2) 

The four fermion interaction can be conveniently studied jlj by introducing a a field in 
the generating functional which satisfies the classical equation of motion a = gipip, and 
transforming as a — > —a under chiral symmetry. The Lagrangian becomes 

£ = tpid^j - \j 2 + gotyty (2.3) 

and by integrating over the fermion field ip one can study the effective action as a func- 
tional of a. 

a) T = 0,/i = 

At zero temperature T and chemical potential /x the calculation indicate the vacuum 
expectation value (a) ^ 0, and therefore the discrete chiral symmetry is spontaneously 
broken. This happens for any value of the coupling constant. There is no contradiction 
with the no-go theorems 0@. 
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b) T ^ 0, /i = 



In this case the time dimension is bounded by (3 = 1/T. The thermodynamical limit 
can be taken only on the space dimension (unless (3 — > oo which goes back to case a)). 
Thus, as far as the occurrence of spontaneous symmetry breaking is concerned, the model 
behaves as one dimensional and the discrete chiral symmetry is immediately restored at 
T c = + . This is a manifestation of the Mermin- Wagner theorem 0. The restoration of 
chiral symmetry at finite temperature is driven by the presence of kinks and antikinks. 
These are non constant field configurations connecting the two degenerate minima, whose 
number grows with the volume IJ. This happens for any large, finite N. If however the 
limit N — > oo is taken before the thermodynamical limit, these configurations are sup- 
pressed, and one is left within the mean field theory, where the model exhibits a second 
order phase transition at a critical temperature T c ^ 0. 

c) T^0,/i^0 

By considering the model in the N —> oo limit, one can derive analytically the phase 
diagram in the plane of chemical potential and temperature (//, T). It turns out that 
there exists a tricritical point separating second order phase transitions from first order 
ones |§. This result has been criticized by a lattice study |10| where the authors find 



that at any /i / the phase transition is first order. They claim that this is due to the 
formation of kink-antikink configurations which are now not suppressed and trigger the 
phase transition in this case. 

In our study of the discrete symmetry model at finite temperature and density a bare 
mass is included, always taken different from zero. This ensures the elimination of kink- 
antikink configurations, which are suppressed in the thermodynamical limit (even in 
presence of a finite chemical potential Q). 

The model has been considered in the iV — > oo limit and we believe that a 1/N expansion 
would not destroy the qualitative picture provided the mass is always kept different from 
zero. 



Continuous symmetry 

We summarize now the main results for the model with continuous symmetry. One starts 
from the Lagrangian 

x One can convince oneself about this fact by looking in ref. || at the procedure leading to the 
expression for the number of kinks, which is independent of the chemical potential. 



C = $i&/> + l -g 2 (^) 2 - l -g 2 (^) 2 (2.4) 
which is invariant under the continuous chiral transformation 

if; -> e i9 ^ 5 ip (2.5) 

In this case one introduces in the generating functional, besides the scalar field a also a 
pseudoscalar field satisfying the classical equation of motion 7r = igip^ip [[[[]. The a and 
7r fields transform under the continuous chiral transformation as 

(a\ ( cos 26 sin 26 



a 



,7r / \ — sin 2$ cos 26 J \7r. 



(2.6) 



Integrating over the fermion fields one gets the effective action as a functional of the two 
fields a and tt. 

a) T = // = 

Due to chiral invariance, the effective potential is a functional of p 2 = a 2 + 7r 2 only. 
Thus one easily realizes that in the mean field approximation its analytical expression is 
equivalent to that of the discrete symmetry model. Then chiral symmetry appears to be 
spontaneously broken. This is also the result at the leading order in 1/N since there is 
no kinetic term for the collective fields. Anyway, rigorously speaking, continuous chiral 
symmetry cannot be broken since there are no Goldstone bosons in two dimensions ||. 
This shows up at the next-to-leading order in 1/N, due to IR divergencies in the ir corre- 
lation function. Nonetheless, as shown in refs. ]7| |12[| , by choosing a polar representation 
for the physical fields, longitudinal and transverse 

a + m = p e id (2.7) 

one finds that the system still allows for an "almost long range order", since the two-point 
correlation function decreases with a power law for large distances, the exponent being 
1/N [00- in deed this behaviour, which was originally found in the low-temperature 
phase of the continuous symmetry X-Y spin model by Kosterlitz and Thouless |I3|, is 



simply dictated by dimensional arguments. Actually it comes out from the fact that one 
is considering the transverse fluctuations of the order parameter, which in the IR are 
dominated by the free propagator of the n field, which is ~ 1/k 2 |l]][12|. Thus the Fourier 
transform in two dimensions diverges logarithmically for large distances (the same holds 
for the 6 propagator). Anyway, by using the physical representation (|2.7|) this behaviour 



has to be exponentiated. As the hypotheses are general, the Kosterlitz-Thouless phase ap- 



pears simply as the "two dimensional version of long range order" in a class of models [14 



A similar phenomenon occurs in 2+1 dimensions for massless bosonic theories at finite 
temperature, as it appears from the behaviour of the free two point correlation function, 
(see Appendix D). 



b) T^O 



Again, for large but finite N, at T 7^ 0, one can easily find that the system behaves as in 
one dimension and that even this "almost long range order" disappears (see Appendix D). 
Anyway, even in this case, the iV — > 00 limit can be formally taken in a way to eliminate 
any space-time dependence of the correlation function. Thus one is left within the mean 
field description, which gives the same critical temperature as in the discrete symmetry 
model. This solution seems to contradict the Mermin- Wagner- Coleman theorems. We 
are going to comment on this point in the following. 



It is evident that whenever non constant field configurations change the picture provided 
by the mean field approach, the 1/N expansion fails. In fact the fluctuations can be 
suppressed only if N is taken strictly infinite, rather than approximating with its leading 
term a 1/JV expansion. Thus, if N is large but finite, the final result is simply that the 
symmetry can be broken spontaneously only if T> > 1 (discrete symmetry) or T> > 2 
(continuous symmetry), T> being the number of dimensions for which one takes the ther- 
modynamical limit. It is clear that going from a zero temperature to a finite temperature 
field theory, T> decreases of one unity (this will be further evidenced in the following and 
in Appendix D). 

To specify something more about the infinite N, let us recall some qualitative arguments 



which lead to the Mermin- Wagner- Coleman results [15[ for a continuous symmetry. A 
criterion to establish the possibility of spontaneous symmetry breaking is to consider the 
ratio 

Z(a) ^ Ar 



Z(0) 



e- m {21 



(r is the euclidean action) between the partition function infinitesimally rotated after 
an operation of the symmetry group, and the unrotated one. To allow for spontaneous 
symmetry breaking, this ratio has to vanish in the thermodynamical limit jnj, and thus 



Ar = T(a) — r(0) has to diverge in the thermodynamical limit. The opposite (i.e. a 
finite Ar) would imply that, being the system initially at one minimum, it would have a 
non zero probability to make a transition to another degenerate minimum. This would 
ensure the order parameter to be zero Actually, only the kinetic term contributes 



to Ar, which at zero temperature implies, in D dimensions 

Ar ~ J(d ll a) 2 d D x ~ L D ~ 2 (2.9) 
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which shows that if D < 2 no spontaneous symmetry breaking occurs. At finite temper- 
ature the time dimension is bounded and the difference in the actions becomes 

AI> ~ J dx Q J {d^afdP^x ~ L°- 3 ~ L d - 2 (2.10) 

where d is the number of spatial dimensions. This phenomenon is much the same as 
the dimensional reduction at high temperatures. There, one is considering the high- 
temperature limit with the spatial dimensions bounded, here the large spatial-dimension 
limit at fixed temperature. Thus a dimensional reduction occurs whenever the ratio 
x'/P^oo (i= 1,2, ....(£> - 1)). 

Finally, if the theory has an internal symmetry group O(N), a N factorizes in the previous 
expressions 

Ar ~ Nf3L d - 2 (2.11) 

and one is lead to discuss in addition the implication of a iV — > oo limit. For instance, 
by allowing iV — > oo as L^, to make Ar divergent for L — > oo it would be enough to take 
£ > 2 — d. Thus if d — 1, the necessary condition is £ > 1, namely that N goes to infinity 
faster than L. 

The large distance behaviour of the free propagator for massless bosons in D — d + 1 
dimensions at finite temperature is dominated by Tx(free massless propagator in d di- 
mensions at zero temperature) (see Appendix D) 

r d d k p* k x 

D^x) ~ T / T^-nr; for |x| - oo (2.12) 



(2yr) d k 2 

Thus, as the IR leading term of the action for transverse fluctuations has the general 
form fl2.9|) , the expression ( |2.12|) leads to the same statement of eq. ( |2.10|) . In presence of 



an internal symmetry group 0(N), a factor 1/N appears in ( p.!2j ), and then, by putting 
an IR cutoff ~ 1/L, one finds that in the IR limit the fluctuations go as 

Dp-j^—; for/^0 (2.13) 
which gives the same information as Q2.ll 



The considerations about the iV — > oo limit are of course not completely satisfactory, as 
the limit appears a bit tricky and not physically clear. However, as the 1/N expansion is 



non-analytic ||12j| , to consider the N oo limit looks rather as moving to a new distinct 
model. 



To summarize, several authors || 0] |TTJ have considered the massless Gross- Neveu 
model in view of its possible similarity to 3+1 QCD taken in the chiral limit. Consequently 
they have neglected the role played by fluctuations, which can destroy the order at low 
dimensionality, but should not be so crucial in 3+1. We also believe that the study of the 



s 



model in the mean field approximation may represent a good guide to the thermodynamics 
of chiral symmetry restoration in 3+1 QCD models. However, both for the discrete and 
continuous symmetry model, we consider important to add a non vanishing bare mass 
term. This allows to escape the Mermin-Wagner-Coleman theorem, and besides it should 
represent a more realistic simulation of the QCD phase transition. We have already 
studied the behaviour of the condensate at finite temperature and density in a massive 
QCD-like theory in the mean field approach |4[], with results qualitatively equivalent to 
those presented here. In this context it is interesting to study further properties of the 
Gross-Neveu model, where the technical complexity is largely reduced. 
We stress again that, differently from QCD, where the zero mass limit can be safely taken, 
in the present case the mass has to be kept finite. In the continuous symmetry model, 
for instance, the illness of the zero mass limit is evident by the drastic change in the long 
distance behaviour of the free correlation function, which passes from an exponential 
decay to a logarithmic (or linear) divergence. We remark that a four fermion model 
in 2+1 dimensions would a priori not be a good candidate to simulate QCD, as chiral 
symmetry can be properly defined only in an even number of space-time dimensions. Thus 
it is not surprising that the phase diagram for symmetry restoration is very different in 
that case [12 . 

In the following we will refer to the mean field results for the massless case only to better 
clarify the results of the massive case. 



3. The massive Gross-Neveu model - Discrete sym- 
metry 

In this section we examine the Gross-Neveu model when the discrete chiral symmetry is 



explicitly broken by adding a fermion mass term to the Lagrangian (|2.1|) 

C = ^d-M)tP + ^g 2 (^) 2 (3.1) 

As already mentioned for the massless case the model is studied by introducing in the 
Lagrangian a a field satisfying the classical equation of motion a = gipip. In presence 
of a mass term it is more convenient to redefine the a field by shifting it by a constant 
a — > o + M/g. Then the Lagrangian (apart from constant terms) reads 

£ = $idi[) - -a 2 + gtjipip - M— (3.2) 
2 9 

Let us first consider the model at vanishing temperature and density. To study the 
problem of spontaneous symmetry breaking (SSB) a suitable tool is the effective potential 
V(a), which has to be minimized with respect to o in order to extract the physically 
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relevant quantities. As anticipated in the previous section, we consider the model in the 
infinite N limit. Following ref . ||16|| , the effective potential is obtained as 

d 2 k 



V(a) 



-a 



iN 



(2tt)* 



\n(g 2 a 2 — k ) + M 



a 
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(3.3) 



The renormalization can be carried out by adding a counterterm 5Za 2 /2 and by imposing 
the condition 



d 2 V/da\-_ 



Then we have 



V(a) 



1 



Xa 2 



2 a2+ An 



In ( % ) 3 

of 



a 

+ M- 
9 



(3.4) 



(3.5) 



The subtraction parameter <jq is arbitrary. A change in <To is equivalent to a change in 
g, to scaling the a field, and to rescaling the bare mass M. Thus the effective potential 
must obey the renormalization group equations 



00 



_d_ 

d(7n 



d d d \ 



(3.6) 



If the potential in eq.(|3.5|) is inserted into eq.( |3.6|) we find 

iV/(27r) 



la 



l+Ng 2 /{2ir) 

9 
2 7. 



(3.7) 



The proportionality between /3(g), jaig) and ^u^g) arises from the fact that to order 
1/N there is no renormalization of the wave function for the ip field at the one loop level. 
Therefore the renormalization of a is determined by the term gaipip, that is from the 
renormalization of g. 

Let us recall that for M = the symmetry is spontaneously broken and the fermion 
acquires a dynamical mass 



m = -go = ga e 



1-tt/A 



(3.1 



where a is the absolute minimum of the effective potential. Eliminating the dependence of 
the effective potential on er by using the relation ( |3.8| ) between a and m , the expression 
ijj) can be put in a manifestly invariant form 



Via) = 



In I 

m J 



M 

g 2 



(3.9) 



in 



We now consider the generalization to finite temperature T and chemical potential fi. 
By using standard techniques, one can derive the effective potential 



V(*,T,ri = V (a) + Vp(<r,T,ii) 



N 2 2 



N 

'J 



m J 



M 

+ B +ga— 



dk 



ln (l + e ~P( V P + 9 2 * 2 + + ( M _> 



(3.10) 



As the pressure equals the negative of the effective potential evaluated at the physical 
points, the equation of state for the system is given by 



P 



-V(T,n,&(T,n,M)) 



(3.11) 



To normalize to vanishing pressure (and energy density, in the limit M — > 0) at zero 
temperature and chemical potential, we have inserted in the expression for the effective 
potential (|0| the "bag term" B = Ng 2 a 2 {0, 0, 0)/4vr (<r(0, 0,0) = a = -m /g). Finally, 
it is convenient to introduce dimensionless quantities 



7T 



m 
r = T/mo 



Nml 



V 



a 



i] = n/m 



y = k/m 



and 



V {u) + Vp(uj,r,r)) 



-cu 2 (ln^ 



1) + - + auo 

' 4 



r — U/ 

-r I dy ln(l + e r ' 
lo L 



v) 



which, at the absolute minimum Co, gives the dimensionless pressure 

7T 



V 



p = -v(u) 



(3.12) 
(3.13) 



(3.14) 



(3.15) 



At zero temperature the effective potential ( |3.14j ) becomes 
V{u) r=0 = V (u) + Vp{u,r = 0, V ) 

= -LU 2 (lntu 2 — 1) H h oluj 

4 v ' 4 

+-6{rf - uj 2 ) [-V\/v 2 -uj 2 + uj 2 ln(y^ _ uj 2 + 77) - -^ 2 lnc; 2 ] (3.16) 

In the following we shall consider the phase diagram and the equation of state for 
small values of the symmetry breaking parameter a. 
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Phase diagram and thermodynamics 

At r = rj = the symmetry is explicitly broken, and the absolute minimum of the 
effective potential is on the negative cu-axis. It appears in correspondence of the point 
Cj ~ —(1 + a), and it becomes deeper and deeper for increasing a, as shown in Fig. 1. 

At finite r and r\ the symmetry still remains broken due to the presence of the breaking 
term au. Nevertheless, if the mass parameter is small (a< 1 but different from zero), 
the system at finite temperature and density exhibits a line of transitions. To be more 
specific, let us start by looking at the behaviour of the condensate in the plane of chemical 
potential and temperature (r],r). To describe the results we find convenient to first 
recall the phase diagram which derives from the study of the absolute minima of the 
effective potential (|3.14|) with a = 0, already found in ref. || (see Fig. 2). There is 
a tricritical point P t = (r] t ,r t ) = (0.608,0.318), separating a line of second order phase 
transitions Lu starting from the point (0, r c ) = (0, 0.567) and ending at P t , from a 
first order line Lj starting from P t and ending at the point (i] c , 0) = (\^2/2,0). In the 
present case, with a ^ 0, by moving along vertical lines in the (f],r) plane, we find 
the following behaviour: for chemical potentials higher than a (mass dependent) value 
r) t (a) the condensate undergoes a finite discontinuity at some (mass dependent) critical 
temperature and drops to a value proportional to the bare mass a. Furthermore these 
critical temperatures decrease for increasing chemical potentials. Thus the system still 
allows for "first order lines" Lj(a) starting from points (r) t (a), r t (a)) and ending at points 
(r/ c (a),0). On the other hand, for r/ < r] t (a), the condensate decreases continuously for 
growing temperatures and thus it is no longer possible to extract a line of second order 
critical points. Nevertheless there still exists a sharp interval of temperatures in which 
the condensate decreases steeply to a small value (proportional to a). This simple result 
can be summed up by saying that by taking into account a small bare fermion mass, the 
phase diagram of the discrete chiral Gross-Neveu model at finite temperature and density 
is a perturbation of the naive mean field solution of ref. §|, found in the chiral limit. 
Although the last one can suffer for the presence of non constant field configurations 
which dominate in the thermodynamical limit, the addition of a bare mass is sufficient 
to eliminate these complications. The 1/N corrections should not destroy the qualitative 
picture in this case. In Fig. 3 we show the comparison between the phase diagram for 
a = 0, for a = 0.01 and for a = 0.1. 

In Fig. 4 we show the condensate behaviour vs. temperature at zero chemical potential for 
a = 0, 0.01, 0.1. We see that for a < 1 the condensate behaviour is still reminiscent of 
a second order phase transition. The same holds for any rj < i] t (a). The low temperature 
behaviour can be derived analytically and it is approximately given by 




(3.17) 



1? 



In Fig. 5 we plot the condensate behaviour vs. chemical potential at zero temperature 
for a = 0, 0.01, 0.1, as derived from the effective potential given by eq. (|3.16|) . In this 
case the critical chemical potential can be evaluated analytically for a < 1. It turns out 
to be 



V2, 



l + 2a) 



(3.18) 



Let us now discuss the equation of state. From eqs. (|3.15|) and ( |3.14|) we have 



p(r,r}) 



UJ 

T 



(l - Mb' 



dy In 



1 

_ 4 

1 + e 



"(V y 2 +& 2 +v)/r 



+ {v -> -v) 



(3.19) 



(we recall that the minimum Cj depends on r, r\ and a). 

The fermion number density and the entropy density are obtained by the derivatives 
of the pressure P with respect to the chemical potential and temperature respectively. 
However it is convenient to define dimensionless number- and entropy-densities 



n 



dp 
dr] 



dy 



1 + e 



[y/v 



iv -> -v) 



(3.20) 



and 



dp 
dr 



y/y' 2 + Cj' 1 + 7] 

[\/y 2 +u 2 +ri]/r 



1 + e 



+ ln(l + 



-\\fy 2 +u 2 +rj\/r 



+ {n^ -v) (3.21) 



which are the true quantities rescaled by a factor TT/(Nm ). Consistently, the dimension- 
less energy density is given by the expression 

dp dp 
—p + r— + T]- 



dr dr] 
—p + rs + rjn 



U) 

T 



(lnt^ 



1 

+ - + acu + 
4 



o i _|_ e [y/y 2 +u 2 +v]/r 



(77 — > —77) (3.22) 



As we have previously done for the phase diagram, it is useful first to refer to the results 
deriving from ( |3.14 ) with a = 0. In this case the exact values of p, n, s, e can be derived 
analytically at r = 77 = and for r > r c (r]) (where r c (rf) is a point of the curves Ljj(a = 0) 
or Lj[a = 0)). In fact, for a = and r = rj = 0, from ( 3.14|) , one has Cj — 1, and thus 
from eqs.( |3.19|) -( |3.22| ) p = s = n = e = 0, whereas for r > r c (r)), one has Cj — 0, and thus 
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the integrals in eqs.( j3.19Q -( j3.22j ) can be exactly evaluated, giving 

7r 2 r 2 n 2 

P = — + \~ b (3 ' 23) 
7r 2 r 

a = — (3.24) 

n = 7] (3.25) 

e = p + 2b (3.26) 

(where b = 1/4 is the dimensionless bag term) . 

In Fig. 6 we show the curves e/r 2 and p/r 2 versus r, for rj = and a = 0. Notice that 
for r > r c ~ 0.57 the two curves approach asymptotically the constant value K = tt 2 /6. 
In fact their expression in this region is given by 

p 7i 2 rj 2 b 



r 6 2r r 

e n 2 r] 2 b 

r 2 6 2r 2 r 



and thus for rj — 0, the two curves in this " Stefan-Boltzmann regime", are symmetrical 
with respect to the dashed line depicted in Fig. 6. Notice also that since both e and p 
are exponentially vanishing for r — ► 0, their ratio to r 2 is vanishing too. 



In Fig. 7 there are the same quantities as in Fig. 6, for rj = 0.63 > r] t (a = 0) ~ 0.608. 
Notice that in this case the curve of the energy density shows a latent heat as the phase 
transition is first order, and that the two curves are no longer symmetrical with respect 
to the asymptotic value K. 

For a ^ 0, the pressure and energy density at r = rj = are corrected to p ~ a and 
e — —a. This is expected because the energy density has to be proportional to the con- 
densate, which means e = au ~ —a since Co ~ —(1 + a) at r = rj = 0. A consequence 
is that at low temperatures the ratios p/r 2 and e/r 2 are now divergent. Thus, to con- 
struct curves analogous to those in Figs. 6 and 7 (where a = 0) for the massive case, 
we have to subtract the divergence. Actually in Figs. 8 and 9 we show [e + a]/r 2 and 
[p — a]/r 2 vs. r at fixed rj with 77 = < T) t (a) and 77 = 0.65 > T) t (a) respectively, for 
a = 0.01. Notice that the latent heat is slightly reduced with respect to the massless case. 

We now come to alternative representations of the phase diagram, different from 
those of Figs. 2 and 3. We consider the planes (l/n,p) and (rt, e), of inverse fermion 
density-pressure and of fermion density-energy density, which are quantities of more direct 
physical interpretation than r and rj. 

Let us start from the isotherms in the (l/n,p) plane for the massless model, in Fig. 
10. These van der Waals type curves are readily obtained in the following way: the 
critical line Lji is the mapping of the corresponding curve in the (77, r) plane (see Fig. 2), 
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whereas the line Lj has split into two parts, corresponding to the value as of the density 
n at the two degenerate minima at the critical chemical potential at a given temperature 
for r < r t . The region inside is the coexistence region. Notice that it is not necessary to 
follow the Maxwell procedure to draw the isotherms, as all the informations are naturally 
included in the formalism. In fact the value of the pressure at the boundaries of the 
coexistence region for a given temperature is directly the negative of the value of the 
effective potential at the two degenerate minima at the critical chemical potential. The 
"liquid phase" corresponds to the region where "the symmetry is restored" and thus where 
the absolute minimum is Co = 0, and the "gas phase" corresponds to the region where "the 
symmetry is broken", or Co 7^ 0. The intermediate values of n are of course not accessible 
from the minima of the effective potential as they only give informations concerning the 
pure phases. They are however physically accessible, in analogy, for instance, with the 
phenomenon of the melting of ice at zero temperature. The supercooled and overheated 
parts of the isotherms can be easily obtained by following respectively the histories of 
the minima at the origin and out of the origin, after they cease to be absolute minima. 
The unstable region is obtained by following the maximum in between, since it appears 
together with one of the two minima until it merges with the other. The picture obtained 
in this way obeys the Maxwell construction which is simply derived from the equality of 
the chemical potentials at the edge of the two phases. 

The isotherms in the symmetric phase have a simple analytical expression for a = 0. In 
fact, from eqs. ( |3.23| ) and ( |3.25| ) which is n = rj in the symmetric phase, one has 

ttV 2 1 

p = 1 — - — — — — b (symmetric phase) (3.29) 

6 2(l/n) 2 

The isotherms in the "gas phase" have been obtained numerically. 

The curve 7 in Fig. 10 is the isotherm at zero temperature, and therefore it defines the 
edge of the accessible region. From the point C = (l/n c , 0) = (a/2,0) it coincides with 
the horizontal axis (up to infinity as n = for rj < rj c ). Thus the line 7 has the simple 
expression 



(3.30) 



In Fig. 11 are shown the isotherms in the plane of inverse density and pressure for 
a = 0.01, as derived numerically. The change in the shape of the coexistence region, 
close to the critical point P t , is well understood by classical arguments of a Landau ex- 
pansion (see Appendix A) [|l7j . 

As one can see, the description for a / (a < 1) is strictly related to that for the 
massless case. This is not surprising as the model admits a Landau expansion which does 
not involve problems related to long range fluctuations if a 7^ 0. Thus the system still 
shows up in two different phases below the critical pressure p t . The situation for a 7^ 
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is much the same as for the vapour-liquid transition for water. It is possible to obtain 
Clausius-Clapeyr on-like equations from the equality of the pressure at the edge of the 
two phases (see ref. [0). 



Finally we come to the phase diagram in the plane (n, e), starting from the case a = 
in Fig. 12. The picture clearly indicates the different regions of chiral symmetry breaking 
and restoration, and the region of coexistence, by the mapping of the corresponding lines 
in Fig. 10. Here the only isotherm drawn is the line 7, which is that for zero temperature, 
and thus it defines the edge of the accessible region (as in Fig. 10). It has the simple 
expression 

«"0H£»+'(-£)G+t) (331) 

The second term in the r.h.s. is obtained from eqs. (|3.23j ), ( |3.25| ) and ( |3.26| ) by taking 
into account that n c = r] c = y/2/2 and 6 = 1/4. As far as the first term is concerned, let 
us first recall that, at zero temperature (see also the second line of eq. (|3.22|) ) 

e = —p + r/n (3.32) 

As the phase transition at r = is first order, the previous formula directly gives the 
latent heat in passing from one phase to the other 

disc e = r] c disc n (3.33) 

(namely the coordinates of the point C = (y/2/2, 1/2)). By allowing for the intermediate 
values of n (which correspond to the physical values in the mixed phase) and taking 
into account that p c (r] c ,r = 0) = 0, whereas the "pure phase" values of n are n\ = 0, 
n 2 = V / 2/2, one gets a straight line 

e = —p c + r] c n = 6 {—^- — n\ n (3.34) 



which is the first term in the r.h.s. of eq. ( |3.31|) . 

In Fig. 13 there is the phase diagram in the (n, e) plane for a = 0.01. Also from this 
picture (as in Fig. 11), apart from the obvious disappearing of the line Ljj, the strong 
effect of the addition of a bare mass on the shape of the curve in the vicinity of the point 
P t becomes evident. 



4. Continuous symmetry 

In this section we consider the massive Gross-Neveu model described by the Lagrangian 
C = $(id - M)ip + \g\^f - ~<?WV0 2 (4.1) 
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which for M = is invariant under the continuous chiral transformation ( j2.5|) . With 
M^O the symmetry is broken from the beginning and there is no reason for the appear- 
ance of exact Goldstone bosons or exact massless particles in general. Actually, as far as 
the phase diagram and the thermodynamics is concerned, it is easy to realize that the 
study of the model in the present case leads to the same results shown in the previous 
section. In fact, after having introduced scalar and pseudoscalar fields and integrated 
over the fermion fields, following the same procedure of the previous section, we obtain 
an effective potential of the form 



M 



V(p 2 ,a,T,p) = V^(p 2 )+ga-; 



where V <y0 \p 2 ) has the same form of (glO|) with M 
minimum conditions 

dV_ 
da 



(p 2 = a 2 + TT 2 ) (4.2) 
= and a <-> p By imposing the 



(<T,7r) = (<r,7r) 



o - dv 

d7T 



(cr,7r) = (cr,7r) 



it comes out that 



2a- 



dV (0) 



dp 2 



+ 



M 



(cr,7r) = (cr,7r) 



(4.3) 



(4.4) 



2tt 



dV (o) 



dp' 



(<x,7r) = ((7,7r) 



(4.5) 



and thus tt = 0, and the effective potential reduces to the one obtained in Section 3. 



Pion decay constant 

Our aim now is to analyse the pion mass and decay constant behaviours in T and p in 
the small mass limit. First of all we notice that the axial current 3^ = ipj^ip formally 
satisfies the PCAC relation 

M 

&>J* = 2—Ti (4.6) 

Moreover, for dimensional reasons, n cannot be identical to the canonical pion field. 
Indeed we expect that the physical pion appears in the Lagrangian with a term Cq^ 2 , 
where Co has the dimension of a mass. Therefore we impose that the renormalized prop- 
agator for ip n has the canonical form 

D^x) = <T^(*K(0)> = / ^ )2p2 l _ m2 f V ' X (4-7) 

Comparing this expression with D n (x) = (Ttt(x)tt(O)) , we obtain that Co is the residue 
of D^{p 2 ) at the pole corresponding to the pion mass 

cl=_ ResD n (p 2 )\ p 2 =m 2 (4.8) 

2 And thus gc Q is invariant under the renormalization group. 
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The divergence of the axial current then becomes 
Defining in the usual way the pion decay constant 



we have 



2Mc n 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



In the soft pion limitQ (at T, \i — 0) the pion mass can be extracted in the following way 

d 2 V 
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mi = -ic 2 D-\U) = cl 



dn 2 



d 2 V 



(4.13) 



where c 2 is the residue of D^ip 2 ) on the pole. 

At finite T, /i, defining m 2 (T,n) as the pole of the thermal two-point Green function 
for the ip n field we obtain 

d 2 V c 2 (T,fi)d 2 V 



mi(T,(i) = -ic 2 {T^)[DMV = <f{T,n) 



8n 2 



c 2 a dipl 



(4.14) 



where c 2 (T, /z) is the residue of D n (p 2 ) at the pole. Using the explicit form of the effective 
potential one has 



Therefore, from eq. ( |4.11 ) we have 



~ 2 c 2 (T^f 



(4.15) 



(4.16) 



ml(T,fi)f 2 (T^) = -A 



M c 2 



-go 



(4.17) 



g 2 c 2 (T, fj,y 

Eq.( f4.17|) is just the Adler-Dashen relation at finite temperature and chemical potential. 



By comparing it to the T = /i = case 



M 



m nfn = — 

g 2 



(4.18) 



3 This is the case in which the it propagator can be approximated as ic 2 ,/ {p 2 — m 2 ). 
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we see that ( 4.18Q is modified by the term Cq/c 2 (T, /i) which depends explicitly on the 
temperature and chemical potential (beside the implicit dependence contained in ga). 

At this point we can obtain m n (T, fi) and f n (T, fi) by calculating explicitly the factor 
c(T, n) . If we remember that the inverse propagator is defined in general as 



D~ (p ) =i(l + 5Z)-U(p 2 



(4.19) 



we can find c(T, /x) by evaluating the self-energy for the it field as given by the diagram 
shown in Fig. 14. Carrying on the calculation as done in App.C we finally arrive at 



1 



N 



l + 2a 2 



d 



dq 



+ 



da 2 J 2u a Ve^K+M) + l e ^K-/x) + 1 



(4.20) 



c(T,/i) 2 Ana 2 

where we have defined ui q = \J q 1 + g 2 cr' 2 . 

We notice that eq. dOOj) can also be verified, at least in the chiral limit, by evaluating 
explicitly the axial current matrix element. Indeed, for m n = 0, we obtain 



ip?U = (0\JH\<P*) 



c(T, n) 



VP ii 



(4.21) 



The behaviours of f„ and in T and /i are shown respectively in Figs. 15, 16 and 
Figs. 17, 18. Notice that, if a = could be taken as physically meaningful, f w could be 
used as a physical signal (it is a measurable quantity) of spontaneous symmetry breaking 
(for T < T c ) and of its restoration (for T > T c ). Unfortunately this is not the case, 
because massless bosons in two dimensions are not allowed. Anyway we notice that for 
small values of a we can still distinguish the transition region where f n decreases very 
rapidly. 

As far as the pion mass is concerned we can say that grows up quickly for T > T c , 
showing that in this case one is outside the range of validity of the soft pion hypothesis 
and the behaviour shown in Fig. 16 keeps only a qualitative value. Nevertheless we ex- 
pect the pion mass to be independent of the symmetry breaking parameter (supposed be 
small) once the symmetry is restored, because in such a case the main parameter would 
be the scale of the theory. On the other hand for T < T c , in the case a = 0, the symmetry 
would be spontaneously broken and the pion would be the associated Goldstone boson; 
the mass that it acquires is thus strongly dependent on the symmetry breaking parameter 
a. Anyway a more precise analysis, that we are not interested to carry on in this pa- 
per, could be carried out by considering the full dependence of the pion propagator on p 2 . 
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5. Discussion and Conclusions 



Study of QCD at finite temperature and density is a formidable problem both for 
analytical methods and for numerical simulations. It may then be of interest to perform 
similar studies on simpler models to get an introductory experience in view of the harder 
QCD problems and to have some hints at the type of phenomena that may be present. 
The most easily treatable problems are in dimensions lower than 4, particularly in 2 
dimensions, one of space one of time. Unfortunately the lowest dimensionality brings out 
peculiar features in itself, which have been known since some time, and are particularly 
expressed by the contents of the so-called Mermin- Wagner Coleman theorem. 

A number of analogies, among them asymptotic freedom, suggest the Gross-Neveu 
model in two dimensions as an interesting candidate in order to carry out orientation 
studies of what might happen to the physical QCD problem at finite temperatures and 
densities. In its original massless form the model falls into the theoretical problems 
caused by its low dimensionality, with direct consequences on symmetry breaking and 
phase structure, on the validity of an expansion in the inverse of the number of flavours 
(1/N expansion), and on the reliability of a mean field approximation. 

The simplest massless Gross-Neveu model has a discrete 75 symmetry, spontaneously 
broken at zero temperature and zero chemical potential, without contradiction of the 
no-go theorems for one-space and one-time dimensions. One expects however symmetry 
restoration as soon as T 7^ 0, for any finite value of N, through a kink-antikink formation 
mechanism. One can consider using a mean field approximation only for the strict N — > 00 
limit, where kinks become suppressed and a nominal phase transition takes place at finite 
T. 

When the chemical potential \i is also non vanishing, the phase diagram for iV — > 00 
shows a structure involving a tricritical point. Our main interest is the massive case, 
for which chiral symmetry is never valid and strictly the no-go theorems do not apply. 
The expected suppression of kink-antikink configurations will restore the validity of the 
large N limit, and one finds that some first order phase transition survives in the T — \i 
plane in spite of the explicit breaking of chiral symmetry, with a critical line apparently 
persisting also for large masses. 

The phase diagram for the massive Gross-Neveu model has been shown in Fig. 3, 
for different values of the mass parameter. Beyond some (mass-dependent) value of the 
chemical potential \i the fermion condensate shows a discontinuity at a critical value of 
T, which becomes smaller for higher fi. One has thus lines of first order transitions. 

There are no second order transition lines for finite mass values, but the condensate 
variation appears as very steep for small masses. 
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Corresponding to the first order transitions one finds a latent heat visible in the 
diagram for the energy density. The phase diagram can, perhaps more physically be 
regarded on the plane (inverse fermion density)-pressure, or alternatively on the plane 
(fermion density)- (energy density). The isotherms in the former plane have the aspect of 
the popular Van der Waals curves of water, showing vapour-liquid transition. These Van 
der Waals curves are here directly provided by the formalism, without having to apply 
Maxwell construction. The end point of the region of phase coexistence is a critical point. 
The massive model admits a Landau expansion, as there are no problems arising from 
possible long range fluctuations. From the equality of the pressure at the borders one 
obtains Clausius-Clapeyron like relations. 

A massless Gross-Neveu model exists also with continuous chiral symmetry. For 
vanishing T and fi uncritical use of the mean field approximation would lead analytical 
results translating those of the corresponding discrete symmetry model. One knows 
however that there cannot be Goldstones in two dimensions (Coleman theorem). This 
is reflected by the appearance of infrared divergent pion correlations at the immediately 
non leading order in the 1/N expansion. An almost long range order, corresponding to 
power law correlation decrease, and reminiscent of the Kosterlitz-Thouless behaviour of 
the X — Y spin model, is however expected, also on simple dimensional arguments. The 
" almost long range order " does not survives at any finite T, even for large N. Only if 
N — > oo is intended as constructive definition of a distinct new model one obtains the 
mean field description. Such constructed model, where fluctuations have automatically 
been suppressed, may not be a consistent approximation to a massless Lagrangian model. 

The N oo properties can be analyzed by a comparison of the partition functions 
after and before symmetry rotation, illustrating the kind of dimensional reduction occur- 
ring in a thermal system, and evidencing at the same time the role played by the infinite 
N. One can also discuss the large N limit of the fluctuations in the infrared, through 
study of the massless thermal propagator. 

From the point of view of constructing some treatable approximation to physical 
QCD, suppressing fluctuations may bring closer to the real situation. Independently of 
this consideration our interest has been centered on massive models, which strictly escape 
the no-go theorems of two dimensions. QCD itself has massive quarks, and the general 
warning when passing to some massless limit is that this is certainly more dangerous in 
any comparable two dimensional model than in QCD itself. Nevertheless we get from our 
exercise the conviction that an overall physical content is included in the N oo limit. 
We had previously studied the chiral phase transition in QCD at finite temperature and 
density in a mean field composite model. The analogy of those QCD results to those 
provided by the N — > oo study of the massive Gross-Neveu model is quite remarkable 
and tends to support the impression of a possibly realistic phase space structure. 

The model with continuous chiral symmetry, broken by the mass term, can easily 
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be led back to the analogous discrete symmetry model to qualitatively discuss its phase 
diagram and its thermodynamics at finite T and fi. A main point of the continuous 
symmetry study has been to discuss the properties of the would-be Goldstone, which 
however here can never actualize itself as a Goldstone because of the low dimensionality 
Nevertheless one can again build up for finite mass a kind of Adler-Dashen relation at 
finite T and fi, and derive the behaviours of the " pion decay constant " and of " pion 
mass ". One cannot, of course, use the pion-constant in this case to retrace the symmetry 
breaking, however we find that traces of the transition remains still visible in the finite 
mass model. In conclusion we can only insist on the overall impression of the usefulness of 
the Gross-Neveu model, in spite of its low dimensionality, to illustrate possible behaviours 
relevant to chiral phases of QCD at finite T and \x. We have carried out here a detailed 
study of such behaviours. 
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A. Expansion of the effective potential for small a 

We show some useful formulas to derive the results of ref. || and to approximate our 
results when the mass parameter a is small. We remember that, whilst the effective 
potential we evaluate is a function which can be safely approximated for small a, to set 
a = is a different matter. In fact for the latter case one should take into account that 
the kinetic term of the effective action produces an ill-defined propagator. Thus, every 
result where a is "physically" taken as vanishing, serves only to recover some known 
result or to possibly suggest what would be expected in a 3+1 dimensional model. 
At finite r and rj the effective potential is given by 



+00 



V(w) 



7T 



2A 



uj 2 (l+5Z)-rY J / dfcla[((2n+l)7rr-i7/) + to 2 + k 2 



+ auj 



(A.l) 



We notice that at any finite r (and/or 77), this expression is analytic in the limit uo — > 
(this is not true at r = 77 = 0, due to a logarithmic divergence in the n = term). Thus 
we can perform a Landau expansion to reproduce the right structure of the effective 
potential near to the origin. We write 
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V(cu 2 ) ~ -au 2 + -bu 4 + -co; 6 + au 
2 4 6 

where the expansion coefficients are given by 



0) 



a = 2 



duo* 



b = 2 



LU = 



d 2 v^ 



c = 



ui=0 



uj=0 



(A.2) 



(A.3) 
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By carrying out explicitly the renormalization, from ( 3.14 ), we obtain 



+00 



1 + 



dk 



A" 



+ 



+ 



1 



- 1 



(A.4) 



.[jfe2+l]f k W k+r l)/ r + 1 e (k-v)/r + l 

This expression can be integrated, as shown in App.B. We have 

a = lnvrr + Re ~^(^~~) ( A ' 5 ) 

Following this procedure also for b and c would lead to very complicated expressions. 
Anyway, by noticing that ultraviolet divergencies are present only in the cj 2 -term, we can 
evaluate the coefficients b and c directly from eq. ([A.l ). Carrying out the integration in 
k, after having derived once, and then deriving again, we find 



1 1 



+00 



1 



1 1 



7r 2 r 2 



Re V 

^ [n + 1/2 - ir)/27rr] 3 8 vr 2 r 2 



ReH 3 



IT] 

2nr 



c- 3 



4 \27rr 



+00 



1 



Re V 

^ [n + 1/2 - z77/27nf 



irr 



(A.6) 



(A.7) 



Using now the relations 
the following results 



and (|A.6|), and referring to Fig. 2 , we can easily obtain 



The equation for the line AP t D (see Treml in ref. [nj)> solution of a = 



r = — exp 

7T 



Re 



\rir 



if, 



IT] 

2rir 



(A.i 



The coordinates of the tricritical point P t , determined by the solutions of the equa- 



tions a = 0, 6 = 



(A.9) 



r t = 0.318 , rj t = 0.608 . 

For 1] = the effective potential becomes (e = u/r) 

^ +00 
V{uJ,r,rj = 0) = -LU 2 (lnuj 2 - 1) + otu - 2r 2 J dk\n(l + e" 



Expanding it in e 2 (see ref. |16|) we arrive at 

V(u) ~ l(\mrr - 7 )cj 2 + \/ oj 4 + auo + . . . 

2 4 871^ H 

where 7 = 0.577. . . is the Euler constant. In this case it is sufficient to terminate 
the expansion at the order u A , because the fact that b is positive assures that the 
potential is bounded from below. For a = the equation a = gives us the critical 
temperature 

m 



(A.10) 



(A.11) 



0.567m, 



7T 







(A.12) 



4 By inserting such results in the expression for c we can verify that c(rt, st) > 0, assuring in this way 
the validity of the previous result. 
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It is evident that the effective potential (|A.1|) , since it admits the expansion ( A. 2 ) with IR 



finite coefficients, leads to the standard mean field critical exponents for the description 
of the system in the vicinity of the line Ljj (up to the tricritical point included). 

We notice that the equations that determine the tricritical point can be obtained by 
another method, based on the analysis of the phase diagram structure, specifically on the 
fact that the line P t E is the envelope of the solutions outside the origin of the stationarity 
equation for the potential. In other words (see Fig. 2), if we consider such solutions, it 
turns out that in the region AP t DO there is one and only one curve for each point (there 
is just one minimum 5 ), while in PtDE there are two solutions for each point (a maximum 
and a minimum). Therefore the place where such points became coincident defines a 
curve, tangent to the former, that is just the line PtE. 

Such a curve can be parametrized by the value u that makes the potential stationary 
in the interval [0, 1] (by the way we notice that uj is a continuous function of T and [/,). 
In particular if F(r,s,u) is the equation for the stationary points (outside the origin), 
the envelope curve is given by 

' F(r,s,u) = 

(A.13) 

dF(r,s,u) _ 
duj 

and the tricritical point can be found by imposing uj = 

These equations (evaluated in uj = 0) are formally equal to the equations a = 0, 6 = 0. 
In fact if we notice that 

Ouj duj 2 duj d(uj 2 ) 2 

it comes out that the equation dV/duj 2 = determines the stationary points outside 
the origin, while uj d 2 V/d(uj 2 ) 2 = is just the derivative of the former. Moreover, the 
second equation must be valid for every value of the parameter uj in the range [0, 1], and 
then it reduces to d 2 V/d(uj 2 ) 2 = 0. 

These two equations tell us that the envelope curve is the place where the curvature in 
the minimum outside the origin changes (the effective potential becomes like a parabola, 
and such minimum disappears), and in the uj = limit they reduce to a = 0, 6 = 0. 

B. Coefficient a of the Landau expansion 



In this section we show how to carry out the integrations in eq. (|A.4|) and thus obtain 



eq.( [A.5 ). The first equation can be transformed according to the following formula 

1 1 

,e^" A + 1 + cy +x + 1 



o 



+oo / 1 I 

dyfiv) —t-t + -^t-t)= (b.i; 



'Due to the symmetry we are just considering what happens, for examples, for u> < 0. 
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+00 I r\ 

dx(f(x + A) + f{x - A))—— - / dx(f(x - A) + /(A - x))- 
o e x + 1 Jo e x 

x 



+ / dxf(x) 
Jo 
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from which we obtain 

"+00 



a = 1 + 



dy- 



lo [y 2 + a] 3 / 2 
The first two integrations give us 



r+oo | r+oo 

/ dy — h P \ dx 

J\ y Jo 



2x 



x 2 - \ 2 e x + 1 



a = ln UJ + 7o dX x 2 -\ 2 e* + l 
The remaining integral can be transformed according to 

2x 1 „ /-+ 00 , 2x / 1 



/•-t-oo 

7o x 2 — A 2 e x + 1 



r+oo 2x 

P / <ix— - — 
Jo 



x 2 — A 2 \e x — 1 e 



2;r 



27(A) - 4J(2A) 



and 



/(A) 



P / dx^—- -± - 

Jo e x - 1 (x + A)(x - A) 



]_ /-+00 

lim- / g?x- 



-0 2 Jo 



+ 



lV(x + A)(x- X + ie) (x + X)(x - A - it) 



Using now the result f20| 

r+oo 



dx- 



x 



1 



io (x 2 + fj 2 )(e x - 1) 2 
and carrying out the limit, we obtain 



In 



P 



2nJ (3 V2tt 



(Re > 0) 



7T 



27(A) -4/(2A) = In - + Re 
A 

from which is easy to get the wanted result. 



2i\) 



iX\ / iX 

77 / ~^\2^ 



C. Pion self-energy 

From Fig. 14 we have 

n(p 2 ) = N{-g) 2 /(-l)tr 



7 - f 7 j 

p + k — m k — m 



2A 



m 2 — k(k + p) 



'fc ((p + fc) 2 — m 2 ){k 2 — m 2 ) 
By using the usual Feynman parameterization we can write 



II(p 2 ) = 2A / da 



m 2 — q 2 + p 2 a(l — a) 
q [q 2 + p 2 a(l — a) — m 2 } 2 



7F-, 



In the soft pion limit we can write 

Tl(p 2 ) =U + p 2 U 1 + 



(C.3) 



where 

The inverse of the pion propagator residue on the pole is just —illi. This can be evaluated 
using the Poisson summation 

l£/[(2n + l)^ + vi] =E(- 1 ) n / ^f(Qo)e in ^-^ (C.5) 

We have 

iX f dq + °° d f dq Q e ™p(q -in) 

= W +2iX J 2^2 ( ~ 1)n(_1) d^2y 27(g 2 + g 2 + m 2 ) (C ' 6) 

The integral in <ig can be performed using the Cauchy theorem, and the sum turns out 
to be just a geometric series. At the end we arrive at 

TT * iN 9 N d [ d<1 1 ( 1 1 ^ 

1 ~ c 2 (T,/i) ~ 4^ + z d^J 2^2^ q { e ^^) + 1 + e ^-M) + lj ( ' } 

D. Bosonic free propagators 

Zero temperature 

We recall some known results at zero temperature. 

In D euclidean dimensions, the free two-point function in coordinates space 

(27r) w k 2 + m 2 

satisfies the wave equation 

(□ - m 2 )D (x, m 2 ) = - f J^L e~ ik -* = -5 D (x) (D.2) 



where 



(2*. 



h 1 — 1 



2fi 



Lorentz invariance requires that D (x,m 2 ) depends only on r = y/x^x^; therefore ( [D.2| 
becomes 

r d 2 D — 1 d 

- dr 2 r dr 

By using the identity 

1 



m 2 



D (r,m 2 ) = -5 D {x) (D.4) 



k 2 + m 2 

in eq. ( P-l|) , one gets 



o 



dae~ a{k2+m2) (D.5) 



1 KD_ x {mr) 

D ^ m ) = J^^Djl [ r (D-2)/2 m (2-D)/2] ( D ' 6 ) 

where use has been made of the representation for the modified Bessel functions 

|ar<? z | < 7r/2; fie z 2 > 
which behave JRJ, for fixed v and z — > 

i^o(^) ~ - In (^\ - 7 

~ I T(z/) (~) " ; (Re > 0) (D.8) 



2 v ' \2 
where 7 = 0.577... is the Euler's constant. 

These formulas allow to recover the well known behaviours of the free massless propaga- 
tors in D=2+l and D= 3+1 at zero temperature 

Do(r) = ~ 
Air r 



For D — 1 + 1, performing the m — > limit with fixed r, from ( |D.ij| ) one finds that the 
expression for the free massless propagator in two (infinite) dimensions 

D (r) = - — ln(r) - — In (—) - (D.10) 



2tt y ' 2vr \2J 2tt 
is logarithmically divergent for large distances. 

Finite temperature free propagators in D=d+1 dimensions 

By working in the imaginary times formalism, the passage from Minkowski metrics to 
Euclidean metrics comes out naturally. In fact one has 

= L ran; P-u) 
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with 



—j3 <t < (3, and 



thus (P = l/T) 



= (Jfeo, k) = (2mnT, k) 
x M = (x ,x) = (-ir,x) (D.12) 



7 +°° r+oo <7 d k 



+00 /-+oo W^i- p —i2imTT+ik-x 

£> /J (r,x)=T £ / 22TU , U 2 (D.14) 



One wants to evaluate the x behaviour in the massless limit for D = 1 + 1, D = 2 + 1 
and D = 3 + 1, and in particular the large distances behaviour in the massless limit. 
To this end it is sufficient to separate out the n = term in eq.( p,13D after having set 
m = 



£>^(r,x) = T / + — — / j— - d e (D.15) 



(2n) d k 2 2tt 2 T7 (2vr) d ^ n 2 + a 2 



where a = kj {2nT) and y = 2txTt. 
The sum of the series is 



cog (ray) _ ?r c%(?r - y)] 1_ 

^n 2 + a 2 2a s/i(a7r) 2a 2 1 ' j 



The large distances behaviours are determined by the IR leading term of (|D.14j ). It is easy 



to see, by expanding ( p,15| ) for small k, that the contribution of the sum ( |D.15| ) is sub- 
leading with respect to the first term on the r.h.s. in ( p.!4|) . Thus the free massless prop- 



agator in D=d+1 dimensions at finite temperature for large distances behaves simply as 
Tx(free massless propagator in d dimensions at zero temperature). 
Let us now evaluate the full x dependence for d = 1,2, 3. 

D=l+1 

In this case 

D ^ X)=T ^J-oo 2^An 2 n 2 T 2 + k 2 + m 2 ^ 

and it is evident that the zero mass limit is badly divergent. Anyway this expression has 
a well defined behaviour in x as it can be easily seen by the fact that its derivative with 



9X 



respect to x is a convergent integral even for m = 0. 

By using the result (|D.15| ) the massless propagator can be put in the form 



, . 1 f+"o dk „ . chlkU - t) 

Dg(T,x) = / — COS(kx) — 

2vrio k K ' sh ^l) 



(D.18) 



As expected the IR behaviour of this integral is the same as that of the n = term of 

( rnsp . 

Let us now evaluate the x derivative of (|D.17| ) 



dDp(r,x 
dx 



dksin(kx) ch —r 



o 2tt sh(k\) 

sh(2irxT) 
~2 ch(2irxT) - cos(2tttT) 



T 



(D.19) 



Integrating again one gets 



Dg(r,x) = C(T) - — In ch{2nTx) - cos(2itTt) 

4:71 L 



(D.20) 



where C(T) is an integration constant which does not depend on x. It is easy to verify that 
in the T — > limit it satisfies (|D.10| ) in the coordinate dependence in two dimensions at 
zero temperature, whereas at finite T, for |x| — ^ oo the leading behaviour is that expected 
for one dimension 



Tlxl 



(D.21) 



D=2+l 

We must evaluate 

^oo d 2 k e -i2nnTT+ik-x 



D p (r,x) = T £ J_ 



oo (2tt) 2 4vr 2 n 2 T 2 + k 2 + m 2 



rp +oo „ +00 u j /i \ ~i2irnTT 



27T n: rr' 00 jo 47r% 2 T 2 + k 2 + m : 
Separating out the n = term (from here on remember that x = |x|), one gets 



T „ , , T 



CO 



D (t, x) = —K (mx) + - V cos(27rnTT)K (xV47r 2 n 2 T 2 + m 2 ) (D.23) 

27F * n=l 



The sum of the series for m = is not a closed form 

oo 

cos(ny)K (nz) = 

n=l 



1 I 

+ 



2 + y 2 2 



7 " "<t) 

1 



+-y\- 

2 ^1 y 2 2 + ( 2vm _ y )2 2vm 

rvi -« 

7T 



n=i \J z 2 + {2-nn + 



2 27m 



^ > 0,yei2(D.24) 



(7 = 0.577.. is again the Euler's constant) 

Identifying y = 2ttTt, z = 2tcTx, and taking into account that, as already seen in ( D.7 ), 
for m — > 



K (mz) ~ —In 



7 



V 2 / 

it follows that for m — > the propagator is given by the expression 
D p (r,x) ^-ln(ra) + ^-ln(T) 

Z7T Z7T 
1 1 



(D.25) 



47r ^ X 2 + r 2 

T_ ^ r 1 1 

4tt^ ^T 2 x 2 + {n - Tt) 2 n 

T_ ~ r 1 1 

4tt ^/T^T(^+T7j2 n 



(D.26) 



In the second line we see the term which survives for T — > 0, which, as expected, corre- 
sponds to the first term of 



For T ^ 0, the dominant term for x —>■ 00 comes from summation of the series. We can 



use the summation formula 22 



N 



N f( x )dx + lf(a) + lf(N) 



+ Rq{N) 



(D.27) 



where a,q,N, are arbitrary integers with a < N and q > 0, and 



B 



29 



(2g)! 



f (2 q -i)( N }_ f m-v( a ) 



(2g-i), 



(D.28) 



This formula holds if f^ 2q \x) is absolutely integrable over (a,N). The first Bernoulli 
coefficients are 



2! ~ 12' 4! 



1 



720 



(D.29) 



It is easy to see that for x — > oo the terms in the first line of ( p. 26 ) are dominant and 
that ln(m) is the only infinite constant for m — > 0. The final result, giving the large x 
dependence in the massless limit, can be put in the form 

T 

D/3{t,x) = -—ln(m)+g[T,x,T) 
In 

47T 

4ir 



where 



P + r + Jx 2 + (/3 + T y 



(3-T+ Jx 2 + (f3- r) 2 



(D.30) 



and thus, for x — > oo 



g(-r,x,T) 

9(t,x,T) 



per x — > oo 
1 1 



47T a/x 2 + r 2 



T 



forT^O 



D p (t,x) ^ ln ( x ) 

which is, as expected, the behaviour typical of two dimensions. 



(D.31) 



(D.32) 



D=3+l 



In D=3+l the massless propagator is well defined even at finite temperature 
D p (t,x) 



pJ^^J-oo (2vr) 3 4vr 2 n 2 T 2 + k 2 



-oo 
+oo 

E 



A; sin(kx) 



-i2nnTT 



2lT 2 X n t^oo ^° 47T 2 T 2 n 2 + fc 2 



(D.33) 



By separating the n = term, and using again the sum of the series in ( p,15| ), one easily 
arrives to the expression 



D p (r,x) 



T 



sh(2irTx) 



Anx ch(2-nTx) — cos(2ttTt) 



(D.34) 



For T — > it is easy to recover the four dimensional result already given in the second of 
( p.8|) , whereas for T^O and x — > oo the leading behaviour is that of three dimensions 



D p (r, x) 



T 

Atcx 



(D.35) 
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FIGURE CAPTIONS 
Fig. 1: Effective potential at T = /i = for various values of a. 

Fig. 2: Phase diagram for the massless Gross-Neveu model. 

Fig. 3: Phase diagram for the massive Gross-Neveu model. 

Fig. 4: Condensate behaviour in r at 77 = for a = 0, a = 0.01 and a — 0.1. 

Fig. 5: Condensate behaviour in rj at r = for a = 0, a = 0.01 and a — 0.1. 

Fig. 6: Plot of e/r 2 and p/r 2 vs r for 77 = at a = 0. 

Fig. 7: Plot of e/r 2 and p/r 2 vs r for 77 = 0.63 > rj t at ct = 0. 

Fig. 8: Plot of [e + a]/r 2 and [p + a]/r 2 vs r for 77 = at a = 0.01. Remember that 
e(r = 0, 77 = 0) ~ —a for a 1. 

Fig. 9: Plot of [e+o;]/r 2 and [p—a]/r 2 vs r for 77 = 0.65 > r) t (a) at a = 0.01. Remember 
that e(r = 0, 77 = 0) ~ —a for a <ti 1. 

Fig. 10: Isotherms of the massless Gross-Neveu model in the (l/n,p) plane. 

Fig. 11: Isotherms of the massive Gross-Neveu model in the (l/n,p) plane for a = 0.01. 

Fig. 12: Phase diagram for the massless Gross-Neveu model in the (n, e) plane. 

Fig. 13: Phase diagram for the massive Gross-Neveu model in the (n, e) plane for 
a = 0.01. 

Fig. 14: Pion self-energy at one loop order. 

Fig. 15: Pion decay constant vs. temperature for a = 0, a = 0.01 and a = 0.1. 

Fig. 16: Pion mass vs. temperature for a = 0.01 and a = 0.1. 

Fig. 17: Pion decay constant vs. chemical potential for a = 0.01. 

Fig. 18: Pion mass vs. chemical potential for a = 0.01. 
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